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A first principles density functional based linear response theory (the so called Density Functional 
Perturbation theory has been combined separately with two recently developed formalism for 
a systematic study of the lattice dynamics in disordered binary alloys. The two formalisms are the 
Augmented space recursion (ASR) and the Itinerant coherent potential approximation (ICPA)Q- 
The two different theories (DFPT-ASR and DFPT-ICPA) systematically provides a hierarchy of 
improvements upon the earlier single site based theories (like CPA etc.) and includes non-local 
correlations in the disorder configurations. The formalisms explicitly take into account fluctuations 
in masses, force constants and scattering lengths. The combination of DFPT with these formulation 
helps in understanding the actual interplay of force constants in alloys. We illustrate the methods 
by applying to a fee FesoPdso alloy. 
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I. INTRODUCTION 

The last thirty years have seen numerous attempts at 
setting up a quantitatively accurate theory of phonons 
in disordered alloys. One of the earliest successful ap- 
proximations was the coherent potential approximation 
[4J (CPA). This approximation was a considerable im- 
provement on the existing theories and, in examples of 
homogeneous disorder, was shown to yield configuration 
averaged Green functions which maintained lattice trans- 
lational symmetry and the herglotz analytical properties 
essential for physical interpretation Despite its 
success, particularly in the electronic problem, the CPA 
was a single-site, mean-field approximation and could 
deal with only diagonal (or mass, in the case of phonons) 
disorder. The phonon problem is specifically difficult be- 
cause, in it, diagonal and off-diagonal disorders are im- 
possible to separate. Moreover, the sum rule satisfied 
between the diagonal and off-diagonal parts of the force- 
constants leads to environmental disorder. That is, a 
configuration fluctuation at a site affects the diagonal 
part of the dynamical matrix at its neighbours. Conse- 
quently, we do not expect the CPA to give an adequate 
description of the phonon problem. That this was so was 
discussed in several papers _ E3j which indicated large 
discrepancies between the CPA predictions and experi- 
mental results. In the electron problem too, whenever 
there was off-diagonal disorder, as in the case of alloys 
with large size difference between its constituents leading 
to local lattice distortions 0] or environmental disorder 
as in the case of alloys with short-range order |3^|-[l5|. 
the CPA was found to be inadequate. 

The hunt for adequate extensions of the CPA was quite 
rigorous during the seventies and eighties 0]-[2(j- Most 



of these generalizations were valid for very special types 
of off-diagonal disorder, which were mostly unphysical, 
or violated translational symmetry and herglotz proper- 
ties. Eventually, three approaches emerged as the most 
successful. Two of them were based on the augmented 
space theorem of Mookerjee [2lJ : the itinerant coherent- 
potential approximation (ICPA) of Ghosh et al Q and 
the augmented space recursion (ASR) of Saha et al j2^| 
and Alam et al [2| ■ The former was an extension of the 
ideas of Mills and Ratanavaraksa [23| and Kaplan et al 
[24| and the latter combined the augmented space tech- 
nique with the recursion method of Haydock et al |25| . 
The third was a very different and rather striking ap- 
proach developed by Rowlands et al |2fJ and Biava et 
al [23 (the non-local CPA or NL-CPA) using the idea of 
coarse graining in reciprocal space originally proposed by 
Jarell and Krishnamurthy |28| . 

More importantly a first principles ab-initio theory of 
phonons in disordered alloys is still lacking. Such a the- 
ory is needed in order to gain a microscopic understand- 
ing of the interplay of force constants in the complex 
phenomenon of phonon excitations. Our aim in this com- 
munication is twofold : First, we shall discuss the simi- 
larities and differences between the two methods (ICPA 
and ASR) based on the augmented space theorem. We 
shall apply both the techniques to identical models of an 
alloy system, FePd, and discuss the comparison between 
their results. Secondly, we shall estimate the dynamical 
matrices from a first-principles approach to the parent 
ordered alloys and compare the ICPA and ASR results 
with experiment. We shall argue that first-principles es- 
timates of the dynamical matrices on ordered versions 
do not yield quantitativly accurate results (in compari- 
son with experiment) for the disordered alloys. We shall 
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propose that we need to go beyond and estimate the dy- 
namical matrices from a model of embedded atoms in a 
fully disordered background. 

The outline of this paper is as follows. In section (II), 
we shall discuss in brief what the augmented space the- 
orem is and its application to the problem of phonon 
excitations in disordered alloys. In section (III), we shall 
briefly describe the salient features of Itinerant coher- 
ent potential approximation (ICPA) and the Augmented 
space recursion (ASR). Section (IV) shall be devoted for 
the details of first principles calculation of force constants 
in alloys. The numerical results shall be discussed in sec- 
tion (V). Finally conclusions are drawn in section (VI) 



II. THE AUGMENTED SPACE THEOREM 



The augmented space theorem has been discussed in 
great detail elsewhere j2j| . In this section we shall intro- 
duce only those salient features which will be required by 
us to understand the mathematics and notations in our 
subsequent discussions. 

For a homogeneously disordered binary alloy the Green 
function may be written as : 



<G(k,w 2 )> = -j- expjik- {R-R')} «(i?|(mw 2 -$)- 1 |i?')» ! 



N r.,< 



m = {m A n R + m B (1 - n R )} 5 RR >, 

* = n R n R , + *f* fl , (1 - n R )(l - n R/ ) + *£ B fl , {n R (l - n w ) + (1 - n R )n R ,} , 

(1) 



Here, R, R' refer to lattice positions, <C ^S> refers 
to configuration averaging over random variables in the 
problem. We should note here that the Fourier transform 
in the first equation may be taken only after the con- 
figuration averaging is carried out. The mass and force- 
constant matrices are matrices in the mode space and for 
systems with one atom per unit cell they are 3x3. {n R } 
are the random site-occupation variables which take val- 
ues 1 and depending upon whether the site labelled by 
R is occupied by A or _B-type of atom. The atom sitting 
at R can either be of type A {n R = 1) with probability 
x or B (n R = 0) with probability y. The augmented 
space formalism (ASF) now introduces the space of con- 
figurations of the set of binary random variables {n R } : 

In the absence of short-ranged order, each random vari- 
able n R has associated with it an operator Mr whose 
spectral density is its probability density : 

p{n R ) = x8{n R - 1) + yS(n R ) 

= -- lim Im(tH I ((n H + iS)I - Mr)- 1 \ U), 

(2) 

where x,y are concentrations of the components A and 
B, Mr is an operator whose eigenvalues 1, correspond 
to the observed values of n R and whose corresponding 
eigenvectors {|1r), |0r)} span a configuration space ip R 
of rank 2. We may change the basis to {| |r), | |r)} 

| U) = {v^|lfl)+W|0 fl )}, 

\ U) = {Vv\1r) - V^\0r)} ■ 



and in the new basis the operator M# corresponding to 
n R is : 

= xV R + yV R + y/xyT^, 

where V\ = |T*><T*|, K = \i R )(i R \ and = 
ltii)(liil + |lij)(Tfil are the projection and transfer oper- 
ators in the configuration space ip R spanned by the two 
basis vectors. 

The full configuration space — Yi R IS then 
spanned by vectors of the form | Ttitl • • ■)• These config- 
urations may be labelled by the sequence of sites {C} at 
which we have a |. For example, for the state just quoted 
{C} = | {3, 5, . . .}). This sequence is called the cardinality 
sequence. If we define the configuration | f | . . . | . . .) as 
the reference configuration, then the cardinality sequence 
of the reference configuration is the null sequence {0}. 

In the full augmented space the operator corresponding 
to n R ; is : 

M R = J®...® Mr® ...<g> J®... S 
The augmented space theorem |2]| then yields : 

«G(k,w 2 )»= (k®{0}|(mw 2 -$f V^W), (3) 

where the augmented k- space basis |k ® {0} ) has the 
form 

(1/VN) exp(-ik -R)\R® {0}). 

R 
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The augmented space operators m and 3? are constructed 
from the original random operators by replacing each 
random variable tir by the operators Mjj. It is an oper- 
ator in the augmented space S = Ti®^. The theorem 
maps a disordered operator described in a Hilbert space 
7i onto an ordered one in an enlarged space S, where this 
space is constructed as the outer product of the space H. 



J 



and configuration space \& of the random variables. The 
configuration space VE* is of rank 2 N if there are N sites 
in the system. Another way of looking at the augmented 
space operators is to note that they are collection of all 
possible operators for all possible configurations of the 
system. 

The augmented space operators may be written as Q : 



off 



< m > i + (m A - m B ) ^ V R ® \(y - x) T l R + 

R 

E E ® [ « » x + *&k {(«-*) + H') + V*y + 4 ; ) } 



(2) 
Hi? 



= E^ 



« »t- e ((y - *)< + v^y^ 1 ) E *rk ((f - ^ 



R'^R. 



* = 



- E { - ® ^ + ^ ® T^) + (j, - x) 2 ^ ® + xyT^ ® 7# } 

I 



(4) 



where, 



(i) _ 



RR' 



V RR' 



2*^ 



The sum rule which connects the diagonal and off- 
diagonal parts of the force-constant matrices has been 
incorporated into the formulation. 



III. THE ITINERANT CPA AND AUGMENTED 
SPACE RECURSION 



Ai A' 
A't A 2 



The partition or downfolding theorem then allows us 
to invert this operator in the subspace spanned by the 
average configurations alone. By the augmented space 
theorem this is the configuration average. If we define 
the operator K as (mw 2 — f>), then using the above par- 
tition : Ki = («m»w 2 -«#>) . The downfolding 
theorem and augmented space theorem together give us : 



The augmented space theorem described in the pre- 
vious section is an exact statement. It is a clever book 
keeping technique to include the effects of disorder fluctu- 
ations in the model of phonons in our random alloy. How- 
ever, it is not an algorithm for the approximate calcula- 
tions of spectral and other physical properties of phonons 
in disordered alloys. For that we have to turn to ei- 
ther mean-field approximations like the CPA and ICPA 
or alternatively to the ASR. The coherent potential like 
mean-field approximations begin with a partition of the 
augmented space into a part which is spanned by the ref- 
erence or null cardinality state |{0}) which we shall call 
the average configuration state and the remaining part 
* - IWXWI spanned by fluctuation states : {|{C})}. 
With this partition, any operator can be written in a 
block representation : 



<G(cj 2 )> = 

F = 
£ = 



(Ki -K't F K')" Pl , 
(G^( w 2 )-£(c 2 ))-^ 
is the itinerator 



K 



-Pi 



(5) 

K't F K' is the self-energy (6) 



Here A _Pl and A~ P2 refer to the inverses of the operator 
A in the subspaces labelled by 1 and 2. This is exactly 
the partitioning idea introduced by Srivastava et al [30|. 
Ghosh et al Q next confined themselves to single fluc- 
tuation states of the type \{R}} and went ahead to self- 
consistently evaluate the self-energy in this approxima- 
tion. Adopting their notation ({R}\ A\{R'}) = A^ (R '\ 
they used translational symmetry in augmented space 
and approximated the self-energy and itinerator F 
within the single fluctuation states : 
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s = ^K' t(jR) F iKKR,) K' iR,) 



RR> 



F(«)(«') = G^ R) 



>RR' 



y^yiRKR") p(R' 



)(R') 



(7) 



(8) 



In going from the equation (5) to (7) all contributions 
to the self-energy of configuration states with more than 
one fluctuations in more than one site have been ne- 
glected. Similarly in going from equation (6) to (8), ma- 
trix elements of the itinerator F between configuration 
states with more than one fluctuation present at a time, 
which corresponds to coherent scattering from more than 
one site have been neglected and such states do not con- 
tribute to F and hence to the self-energy £ within this 
approximation. The second equation is a Dyson equation 
within the subspace spanned by only single fluctuation 
states. Self-consistency is achieved through : 



= K ' iR,) F( fl 'H R ") K' {R "l 

R'R"^R 

The above argument shows that unlike the usual CPA 
where only a single fluctuation at a site is considered, 
multiple fluctuations coming from multiple-scattering is 
present in the itinerator F and therefore contribute to 
the self-energy S. However, the approximation described 
above means that correlated fluctuations between more 
than one site are present in neither the itinerator nor the 
self-energy. This is the main approximation involved in 
the ICPA. 

We note that the ICPA is a self-consistent mean-field 
approximation for the self-energy which relates the con- 
figuration averaged Green function to the virtual crystal 
one. It is an approximation which maintains both the 
translational symmetry of the configuration average and 
its herglotz analytic properties. The ASR is an alterna- 



tive technique for doing the same thing, namely obtain- 
ing an approximation to the self-energy maintaining the 
necessary properties of the exact case. 

The recursion method addresses inversions of infinite 
matrices |3^ . The average Green function in the aug- 
mented space formalism can be written as [2J : 



«G»={l|(^' 2 i-K)- 1 |l} 



, Here 

J 
K 

11} 



Where /i 1 / 2 =< to -1 > 



/^ V V /2 ) |{0}> + (^ 1/ V /2 )|{i?}>, 

<C mT 1 !" 1 and 



/'2 



1/2 



xy{(m 



-1/2 



; Mi 

-1/2, 



Once a sparse representation of an operator in a 
Hilbert space, K, is known in a countable basis, the re- 
cursion method obtains an alternative basis in which the 
operator becomes tridiagonal. This basis and the rep- 
resentations of the operator in it are found recursively 
through a three-term recurrence relation : 

\u n +i} = K|u„} - a n (k)\u n } - /3 2 (k)|w„_i}. (9) 

with the initial choice \u\} = |k(g> {1}) or and 1 = 1. 
The recursion coefficients a n and f3 n are real and are 
obtained by imposing the ortho-normalizability condition 
of the new basis set as : 



{u n \K\u n } 



{u n \K\u n+ i\ 



_ L "'"" r '" ; /3 2 +1 (k) — ill 
{u n \u n \ {u n \u„} 

and also {wm|u n } = for m^n,n±l. 

Now, we use the augmented space theorem and re- 
peated applications of the downfolding theorem on the 
tri-diagonal representation gives : 



< G(k,Lu' 2 ) > 



J 2 -ai(k) - 



/3 2 2 (k) 



cj ,2 - ai (k)-E'(k,cj' 2 )' 



(10) 



J 2 - a 2 (k) 



uj' 2 -a N (k)~T(k,uj' 2 ) 



From the definition of the self-energy given earlier, it 
has been argued by us in an earlier communication that 
the disorder scattering induced lifetimes come entirely 
from the imaginary part of £'(k, u/ 2 ). Here T(k, u/ 2 ) 



is the asymptotic part of the continued fraction. The 
approximation involved has to do with the termination 
of this continued fraction. The coefficients are calcu- 
lated exactly up to a finite number of steps {a n , (3 n } for 
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n < N and the asymptotic part of the continued frac- 
tion is obtained from the initial set of coefficients using 
the idea of Beer and Pettifor terminator 

|H. With this 

terminator, the approximate Green function maintains 
the herglotz properties of the exact result. Haydock and 
co-workers [23 have carried out extensive studies of the 
errors involved and precise estimates are available in the 
literature. Haydock |25| has shown that if we carry out 
recursion exactly up to N steps, the resulting continued 
fraction maintains the first 2N moments of the exact re- 
sult. 

Both the ICPA and the ASR involve approximations 
of the self-energy. We have already discussed that in the 
ICPA contributions of configurations involving correlated 
fluctuations in more than one site to the self-energy are 
ignored in the present case but it is capable of incorporat- 
ing them. If we use the form of Eqn. (4) in the recursion 
equations (9), it is immediately obvious that in the ASR, 
contributions of such correlated fluctuation states to the 
self-energy are present. We had earlier shown that such 
contributions occur first at /3| for diagonal disorder and 
in a2 in case of off-diagonal disorder. Ignoring such con- 
tributions will make all moments greater than or equal 
to eight to be non-exact for digonal and three for off- 
diagonal disorder. The ICPA achieves accuracy through 
self-consistency in the subspace of single fluctuations in 
it's present version , while the ASR achieves accuracy by 
increasing the number of recursions in the full augmented 
space and estimating the terminator to mimic the asymp- 
totic part of the continued fraction as closely as possible. 
The two are very different algorithms. Both can take care 
of off-diagonal disorder and short-ranged order, but in 
those situations where clustering, either chemical or sta- 
tistical, is important, the ASR, which takes into account 
correlated scattering from clusters, should be preferable 
over the single-fluctuation only version of the ICPA. 



IV. FIRST PRINCIPLES CALCULATIONS OF 
FORCE CONSTANTS IN ALLOYS 

As is clear from the above discussions that the crucial 
component in both the ASR and the ICPA is the alloy 
force constants. Due to the random chemical environ- 
ment around each atom in a substitutionally disordered 
alloy, the force constants corresponding to A-A,B-B and 
A-B pairs in a A^Bi-x alloy are different and in no way 
resembles the force constants in a completely ordered en- 
vironment. In order to have significant accuracy in calcu- 
lated phonon properties one should, therefore, have ac- 
curate information on force constants corresponding to 
various pairs of chemical species. The only trustworthy 
source of force constant data is the first-principles cal- 
culation. To this end, we have employed first-principles 
Density functional perturbation theory (DFPT) to ob- 
tain force constants. The details about the DFPT and 
our approach to use it to extract random alloy force con- 
stants is discussed below. 



A. Density functional perturbation theory 

Density functional perturbation theory (DFPT) [l| is 
a density functional theory (DFT) based linear response 
method to obtain the electronic and lattice dynamical 
properties in condensed matter systems. The dynami- 
cal matrix which provides information on lattice dynam- 
ics of the system can be obtained from the ground-state 
electron charge density and it's linear response to a dis- 
tortion of the nuclear geometry |34|. In DFPT, this lin- 
ear response is obtained within the framework of DFT. 
One of the greatest advantages of DFPT-as compared 
to other nonperturbative methods for calculating lattice 
dynamical properties of crystalline solids (such as the 
frozen-phonon or molecular dynamics spectral analysis 
methods) -is that within DFPT the responses to per- 
turbations of different wavelengths are decoupled. This 
feature allows one to calculate phonon frequencies at ar- 
bitrary wavevectors avoiding the use of supercells and 
with a workload that is independent of the phonon wave- 
length. 



B. Random alloy force constants from DFPT 

Since there is no first-principles theory for lattice dy- 
namics in random alloys available we took recourse to 
calculate force constants for ordered structures which can 
suitably mimic the random alloy using DFPT. However, 
for a proper representation for the random alloy, one 
needs to work with a large supercell which prohibits the 
use of DFPT from a practical point of view. The other 
approach would have been to construct a set of ordered 
structures having the same composition of the alloy un- 
der investigation, run first-principles calculation on each 
of them and average the data appropriately. As a first ap- 
proximation to this approach, here we have done DFPT 
calculations on a single ordered structure and used the 
resultant force constants as approximate random alloy 
force constants as inputs to the ICPA and the ASR. The 
alloy chosen is FePd. The reason for choosing the FePd 
system is twofold: first, the ICPA and the ASR were 
applied only for the NiPt and NiPd alloys where the con- 
stituents of the alloys have face-centred cubic structure 
in their elemental phases. In case of FePd, although the 
alloy in the disordered phase is face centred cubic but 
Fe is body centred cubic in it's elemental phase. It was 
therefore interesting to test the suitability of both the 
approximations in case of such a system where one of the 
constitutents forming the alloy has a different structure 
than the alloy itself in it's elemental phase, second, in- 
elastic neutron scattering data was available for FesoPdso 
|4l|. It would therefore have been possible to compare 
the ICPA and the ASR results with the experimental 
data directly enabling the understanding of the nature 
of interactions between various pairs of species in the 
random phase. Since the FesoPdso forms a face centred 
cubic (fee) solid solution, we have chosen the prototype 
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tetragonal Llo structure with c/a ratio equal to unity to 
be used for first-principles calculations. 



C. Details of first- principles calculations 

We use DFPT within the local-density-approximation 
(LDA) to compute the force constants for the FePd 
equiatomic composition single ordered structure men- 
tioned above. The experimental lattice constant a = 7.24 
a.u. is used in the calculations. We employ a plane- 
wave pseudopotential implementation of the DFPT with 
Perdew-Zunger parametrization of the LDA |35| as done 
in the Quantum- Espresso package [2B|. Ultrasoft pseu- 
dopotentials with non-linear core correction [2a| are 
used for Fe and Pd. The kinetic energy cut-off is taken to 
be 35 Ry. The Brillouin-zone inte grat ions are carried out 
with Methfessel-Paxton smearing 39] using alOx 7x 7x 
k-point mesh, which corresponds to 120 k-points in the 
irreducible wedge. The value of the smearing parameter 
is 0.1 Ry. These parameters are found to yield phonon 
frequencies converged to within 5 cm -1 . 

Once adequate convergence is achieved for the elec- 
tronic structure, the phonon force constants are obtained 
using the linear response. Within DFPT, the force con- 
stants are conveniently computed in reciprocal space on 
a finite q-point grid and Fourier transformation is em- 
ployed to obtain the real-space force constants. The 
number of unique real-space force constants and their 
accuracy depend upon the density of the q-point grids: 
the closer the q-points are spaced, the accurate the force 
constants are. In this work, the dynamical matrix is com- 
puted on a 6x 6x 4 q-point mesh 40] commensurate with 
the k-point mesh. 



V. RESULTS AND DISCUSSIONS 

In Table I we report the nearest-neighbor force con- 
stants for the artificial ordered structure A-B obtained 
from first-principles as described above. Subsequently, 
we use these force constants as inputs to the ICPA and 
the ASR calculations. Figure QJ shows the correspond- 
ing dispersion curves. The results clearly show that the 
force constants for the artificial ordered structure are not 
adequate to describe the lattice dynamics for the disor- 
dered Fe5oPd5o system. In both the ICPA and the ASR, 
the high frequency phonons are poorly represented for 
all three symmetry directions. On top of that, the high 
frequency branches suffer a split for large q values, a fea- 
ture not observed in the experiments. All these features 
point to the fact that the force constants used in the 
ICPA and the ASR calculations completely fail to cap- 
ture the complexities of the force-constant disorder in a 
random environment. This is quite understandable as we 
have used a crude approximation for the force constants 
in random environment. One single ordered structure , 
in no way, can mimic the randomness in the environment 



TABLE I: Real-space nearest neighbor force constants for 
FesoPdso obtained by DFPT calculations on the artificial or- 
dered structure. The units are dyn cm -1 . 



Pair 


Force constant 


Direction 


Fe-Fe 


-9458 


lxx 


Fe-Pd 


-9458 


lxx 


Pd-Pd 


-28974 


lxx 


Fe-Fe 


-6005 


lxy 


Fe-Pd 


-10755 


lxy 


Pd-Pd 


-30372 


lxy 


Fe-Fe 


1800 


Izz 


Fe-Pd 


-114 


lzz 


Pd-Pd 


3555 


Izz 



around a given chemical species. The fact that the force 
constants obtained on this artificial structure are respon- 
sible for the disagreement with the experiment is corrob- 
orated by the coherent structure factors, particularly at 
high q values as demonstrated in Figure where coher- 
ent structure factors for certain high q values, obtained 
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FIG. 1: Dispersion Curves ( frequency v vs. reduced wave 
vector £ ) for FesoPdso alloy. The upper panel correponds to 
the ASR results however the lower panel to the ICPA results. 
The filled circles are the experimental data EUl . The force 
constants used are given in Table I. 
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FIG. 2: Partial and total structure factors calculated in the 
ICPA for various £ values along the [£, 0, 0] and [£, £, Q direc- 
tions in FesoPdso alloy. The type of mode is labelled along a 
particular symmetry direction. The force constants used are 
given in Table I. The details are given in the text. Descrip- 
tion of various curves 
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by the ICPA, are displayed. The curves clearly show 
that the spurious high frequency peak is due to the Pd- 
Pd pairs and to a smaller extent due to Ni-Pd pairs. As 
is seen from Table-I, the Fe-Fe and Pd-Pd force constants 
differ by about 70 %, thereby representing a situation of 
very strong disorder as is seen in the case of Ni 50 Pt5o 0. 
The splitting of the high frequency branch is a manifesta- 
tion of this strong force-constant disorder, albeit wrong 
in the present case. In the pursuit of the correct set of 
force constants for the system so that the suitability of 
the ICPA and the ASR can be properly tested, we then 
used the force constant data as reported in the experi- 
ment |4l|. However, the force constant data reported in 
the experiment was obtained by fitting the frequencies to 
a Born-Von-Karman force constant model. Also, the fre- 
quencies were obtained from neutron-scattering data on 
ordered L10 structure at 860K, very close to the order- 
disorder transition temperature 950K. The reason behind 
using the experimental force constants obtained from an 
ordered L10 structure for disordered calculations were 
twofold: first, the L10 force constants should be a bet- 
ter approximation for random alloy force constants than 
the artificial cubic structure ones because the L10 struc- 
ture allows structural relaxation and therefore a vari- 
ation in the bond distances between different chemical 
species pairs although in a restricted way. Nevertheless, 
this restricted degree of relaxation could be crucial in 
capturing the nature of forces between various chemi- 
cal species as has been seen in case of NiPt alloys |42j . 
Second, since the L10 data was taken at 860K and the 



FIG. 3: Dispersion Curves ( frequency v vs. reduced wave 
vector £ ) for FesoPdso alloy. The upper panel correponds to 
the ASR results however the lower panel to the ICPA results. 
The filled circles are the experimental data |4l|. The error bar 
in the ASR result basically represents the full widths at half 
maxima (FWHM) at various C, values. The force constants 
used are given in Table II. 

disordered fee data was taken at 1020K, both of them 
lie very close to the order-disorder transition tempera- 
ture. At a first-order order-disorder transformation at 
finite temperatures, the ordered phase is only partially 
ordered and the disordered phase is in equilibrium with, 
has short-range order. Examination of the correlation 
functions has shown that ordered and disordered states 
rather show similar atomic arrangement in the vicinity 
of the order- disorder transformation. It is therefore ex- 
pected that in the present case, the L10 force constants 
at 860K would not change significantly in the disordered 
phase at 1020K. These intuitive arguments are well sup- 
ported by the dispersion curves presented in Figure J3J). 
Both the ICPA and the ASR results agree reasonably 
well with the experimental data. The spurious splitting 
obtained earlier disappears. This disappearance can be 
understood better if we look at the force constants used 
for this calculation. Table-II lists the experimental force 
constants used as inputs for the ICPA and the ASR cal- 
culations. 

In comparison with DFPT values, the Pd-Pd force 
constants are a lot softer and the Fe-Pd force constants 
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TABLE II: Real-space nearest neighbor force constants for 
FesoPdso obtained from experimental data [4l| on L10 struc- 
ture at 860K. The units are dyn cm -1 . 



Pair 


Force constant 


Direction 


Fe-Fe 


-5650 


lxx 


Fe-Pd 


-14050 


lxx 


Pd-Pd 


-19450 


lxx 


Fe-Fe 


-9750 


Ixy 


Fe-Pd 


-16550 


lxy 


Pd-Pd 


-22350 


Ixy 


Fe-Fe 


4100 


Izz 


Fe-Pd 


2500 


lzz 


Pd-Pd 


2900 


Izz 



harden. The fact that the force constants and their be- 
havior is indeed the deciding factor is again exemplified 
by the coherent structure factors for the selected high q 
vectors as shown in Figure 10} . The figures show that 
the single high frequency peak is now mostly because of 
the Fe-Fe and Fe-Pd contributions, rather than the Pd- 
Pd contribution. This points to the fact that the Pd-Pd 
contribution was overestimated and the Fe-Fe and the 
Fe-Pd contributions were grossly underestimated by the 
DFPT calculations on the artificial cubic structure be- 



cause of the lack of relaxation in such structure. This in 
turn can be understood by looking at the bond distances 
between various pairs of species. In the cubic structure, 
the Fe-Fe, Pd-Pd and the Fe-Pd distances were same and 
in the present case was taken to be 5.12 a.u. The L10 
structure at 860 K, on the other hand, had Fe-Fe and 
Pd-Pd distances to be 5.25 a.u. and the Fe-Pd distances 
to be 5.08 a.u. Thus, the Pd atoms, in the artificial 
cubic structures, were made to vibrate in a smaller vol- 
ume and because of the smaller distance between two like 
atoms, the Pd-Pd force constants became harder. Fig- 
ure JSJ) compares the ICPA and the ASR results for the 
phonon densities of states. Both the approximations pro- 
duce identical features. The peaks and the band edges 
have quantitative agreement among themselves and with 
the experimental results 41]. 

Figure © displays the Full width at half maxima 
(FWHM) data associated with the finite lifetimes of 
phonons due to disordered scattering. The upper three 
panels show the widths as a function of wave vector (() 
along the three symmetry directions extracted from the 
ASR method, however the lower three panels show the 
ICPA results. The FWHM is a more sensitive test for 
the underlying approximation than the dispersion curves. 
The ICPA and the ASR have reasonable agreement re- 
garding FWHM's although quantitative agreement is not 
there because of the different nature of the two ap- 



Experimental force constants 




6 8 2 

Frequency V(THz) 

FIG. 4: Partial and total structure factors calculated in the 
ICPA for various £ values along the [£, 0, 0] and [(, £, Q direc- 
tions in FesoPdso alloy. The type of mode is labelled along a 
particular symmetry direction. The force constants used are 
given in Table II. Other details are given in the text. De- 
scription of various curves 




Frequency (THz) 



FIG. 5: Phonon density of states for FesoPdso alloy. The up- 
per and lower panel shows ICPA and ASR results respectively. 
The force constants used are that of experimental paper |41| 
given in Table II. 
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proximations. Unfortunately, the experimental group 
didn't perform a phonon life-time measurement so that 
the FWHM's calculated theoretically could be compared 
with the experiments. 



FIG. 6: Disorder-induced FWHM's vs wave vector (£) for 
FesoPdso alloy. The upper and lower panels show ASR and 
ICPA results respectively. The force constants used are given 
in Table II. 



In Fig. Q, we have shown the total coherent structure 
factors at various £ values along the three symmetry di- 
rections with various modes of vibrations. The upper 
box shows ASR result however the lower ICPA result. In 
both the cases, the different curves for different ( val- 
ues are shifted along the x-axis in order to facilitate vi- 
sion. One can easily notice the difference in the nature 
of curves arising out of two different methodologies. In 
the ASR result, we have found three classes of degener- 
ate modes, these are (1) Tl and T2 modes along [C00] 
directions (2) L and Tl modes along l((0] directions and 
(3) all the three modes along [CCC] directions. However 
in the ICPA case, all the modes along different directions 
are non-degenerate. The difference in the two kinds of 
results are due to the different structure and way of cal- 
culations in the two methodologies. However the ulti- 



mate dispersion curves and FWHM's came out from two 
different methodologies are reasonably well compareable, 
which are the actual vibrational properties to look at in 
any disordered alloy. 



VI. CONCLUSIONS 

This paper has continued the development of the Aug- 
mented space recursion 

2] and the Itinerant coherent po- 
tential approximation (ICPA) 3] for studying the vibra- 
tional properties of disordered metallic alloys. A brief 
description of the two methods combined with a first- 
principles calculation (the so called DFPT) of the dy- 
namical matrices has been reported. The power of these 
approaches has been illustrated by explicit calculations 
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FIG. 7: Total coherent structure factors in different directions with different branches for FesoPdso alloy. The upper and lower 
box shows ASR and ICPA results respectively. In each of the different directions and branches, the various curves indicate the 
total structure factors for various £ values starting from the lowest value to the edge of the Brillouin zone. In the ASR result, 
the Tl and T2 modes are degenerate along [£00] direction, L and Tl modes are degenerate along the [££0] direction, however 
all the three modes are degenerate along the [CCC] directions. In the ICPA result, non of the modes are degenerate. Such 
differences in the two methods are hidden inside the structure of two different calculations. In both the boxes, the different 
curves for different £ values are shifted along the x-axis in order to facilitate vision. 



on the FesoPdso alloy. 

Both the theories are unique and systematic in the 
sense that they produce almost identical results for a par- 
ticular system. Both the theories can explicitly take into 
account the fluctuations in masses, force constants and 
scattering lengths. Wc propose the methods as compu- 
tationally fast and accurate techniques for the study of 
lattice dynamics of disordered alloys. A correct quanti- 
tative trend (compareable to the experimental results) of 
the phonon dispersion and the phonon density of states 
has been predicted by both the methodologies when the 
experimental force constants has been used in the calcu- 
lation. Off-course there is a fairly obvious general com- 
ment to be made with regard to the self consistency of 



the procedure. This is precisely the reason that a first- 
principles estimate of the dynamical matrices on parent 
ordered alloys do not yield quantitatively accurate results 
(in comparison with the experiment) for the disordered 
alloy. We shall propose that we need to go beyond and 
estimate the dynamical matrices from a model of embed- 
ded atoms in a fully disordered back ground. 



Acknowledgment 

One of us (AA) would like to acknowledge the CSIR 
for financial assistance in the form of a Fellowship Grant. 



11 



[1] S. Baroni, S. De Gironcoli, A. Dal Corso and P. Giannozzi 

Rev. Mod. Phys. 73, 515 (2001). 
[2] A Alam and A Mookerjee, Phys. Rev. B 69, 024205 

(2004). 

[3] S. Ghosh, P.L. Leath and M.H. Cohen, Phys. Rev. B 66, 
214206 (2002). 

[4] D.W. Taylor, Phys. Rev. 156, 1017 (1967). 

[5] A complex function f(z) is called herglotz if (i) all its 
singularities lie on the real z-axis (ii) Im f(z) is nega- 
tive everywhere in the plane Im (z) > and is positive 
everywhere in Im (z) < (iii) lim z= E^±ao ^tef(E) = 
0, E is real. 

[6] Miiller-Hartmann. 

[7] N. Kunimoto, Y. Tsunoda and Y. Hirai, Solid State Com- 

mun. 13, 495 (1973). 
[8] T. Kaplan and M. Mostoller, Phys. Rev. B 9 , 353 

(1974). 

[9] W.A. Kamitakahara, Bull. Am. Phys. Soc. 19, 321 
(1974). 

[10] W.A. Kamitakahara and D.W. Taylor, Phys. Rev. B 10 
, 1190 (1974). 

[11] M. Mostoller, T. Kaplan, N. Wakabayashi and R.M. 

Nicklow, Phys. Rev. B 10 , 3144 (1974). 
[12] H.G. Smith and N. Wakabayashi, Bull. Am. Soc. 21, 410 

(1976). 

[13] T. Saha, I. Dasgupta and A. Mookerjee, J. Phys. : Con- 
dens. Matter 8, 1979 (1996). 

[39] A. Mookerjee and R. Prasad, Phys. Rev. B 48 , 17724 
(1993). 

[15] D. Paudyal, T. Saha-Dasgupta and A. Mookerjee, J. 

Phys. : Condens. Matter 16, 2317 (2004). 
[16] H. Shiba, Prog. Theor. Phys. 40, 942 (1968). 
[17] T. Kaplan and M. Mostoller, Phys. Rev. B 9 , 353 

(1974). 

[18] B. Nickel and W.H. Butler, Phys. Rev. Lett. 30, 363 
(1973). 

[19] F. Ducastelle, J. Phys. C 7, 1795 (1974). 
[20] A. Gonis and JW. Garland, Phys. Rev. B 18 , 3999 
(1978). 

[21] A. Mookerjee, J. Phys. : Condens. Matter 6, 1340 (1973). 
[22] T. Saha, I. Dasgupta and A. Mookerjee, Phys. Rev. B 

50 , 13267 (1994). 
[23] R. Mills and P. Ratnavararaksha, Phys. Rev. B 18 , 

5291 (1978). 



[24] T. Kaplan, P. L. Leath, L. J. Gray and H. W. Diehl, 

Phys. Rev. B 10 , 1200 (1980). 
[25] R. Haydock, V. Heine and M.J. Kelly, J. Phys. C 5, 2845 

(1972). 

[26] D.A. Rowlands, J.B. Satunton, B. L. Gyorffy, E. Bruno 
and B. Ginatempo, Phys. Rev. B 72 , 045101 (2005). 

[27] D.A. Biava, S. Ghosh, D.D. Johnson, W.A. Shelton and 
A. Smirnov, Phys. Rev. B 72 , 113105 (2005) 

[28] M. Jarrell and H. R. Krishnamurthy, Phys. Rev. B 63 , 
125102 (2001). 

[29] A. Mookerjee, Electronic Structure of Alloys, Surfaces 

and Clusters edited by A. Mookerjee and D.D. Sarma 

(Taylor Francis, London, 2003). 
[30] V. Kumar, A. Mookerjee, and V. K. Srivastava, J. Phys. 

C 15, 1939 (1982). 
[31] S. Ghosh, N. Das and A. Mookerjee, Int. J. Mod. Phys. 

B 21, 723 (1999). 
[32] R. Haydock Solid State Physics (Academic Press, New 

York, 1980), Vol. 35. 
[33] N. Beer and D.G. Pettifor, in Electronic Structure of 

complex systems, edited by P. Phariseau and W.M. Tam- 

merman (Plenum, new York, 1984), P- 769 
[34] R. Pick, M.H. Cohen and R.M. Martin Phys. Rev. B 1 

, 910 (1970); P.D. De Cicco and F.A. Johnson Proc. R. 

Soc. London, Ser. A 310, 111 (1969). 
[35] J.P Perdew and A. Zunger, Phys. Rev. B 23 , 5048 

(1981). 

[36] S. Baroni, A. Dal Corso, S. de Gironcoli, and P. Gian- 
nozzi, http: / / www . Q uant um-Espresso.org 

[37] D. Vanderbilt, Phys. Rev. B 41 , 7892 (1990). 

[38] S.G. Louie, S. Froyen and M.L. Cohen, Phys. Rev. B 26 
, 1738 (1982). 

[39] M. Methfessel and T.A. Paxton, Phys. Rev. B 40 , 3616 
(1989). 

[40] P. Gianozzi, S. de Gironcoli, P. Pavone and S. Baroni, 

Phys. Rev. B 43 , 7231 (1991). 
[41] T. Mehaddene, E. Kentzinger, B. Hennion, K. Tanaka, 

H. Numakura, A. Marty, V.Parasote, M.C. Cadville, M. 

Zemirli and V. Pierron-Bohnes, Phys. Rev. B 69 , 

024304 (2004) 

[42] S.Ghosh, J.B. Neaton, A. H. Antons, M.H. Cohen and 
P.L. Leath Phys. Rev. B 70 , 024206 (2004) 



